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ABSTRACT 



An analysis of fin-line characteristics is presented. 

The dispersion characteristic and characteristic impedance 
of the structure are determined by the application of the 
spectral domain technique. Numerical results are presented 
and compared with known data for ridged waveguide, dielectric 
slab loaded waveguide and slotline. Excellent agreement has 
been obtained. Design curves are also presented for practi- 
cal millimeter wave fin-line structures covering the 26.5 - 
140 GHz frequency range. 
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INTRODUCTION 



I . 



A . BACKGROUND 

A variety of waveguiding structures have been proposed 
for use at microwave and millimeter wave frequencies [Ref. 1] . 
In recent years, one of these structures, fin-line, has gained 
in importance as a transmission medium in millimeter wave 
circuit construction [Refs. 2-7]. Fin-line has been found 
superior to microstrip at millimeter wavelengths as the 
former provides eased production tolerances, better compati- 
bility with hybrid devices, greater freedom from radiation 
and higher mode propagation, combined with the ability to 
construct simple transitions to conventional rectangular 
waveguide. 

Figure 1 shows a 3-dimensional view of fin-line. This 
structure is interesting in that it may be viewed in various 
ways depending upon the value of W/b, the ratio of slot width 
to waveguide height. For small values of W/b, the structure 
may be appropriately viewed as a slotline with a rectangular 
shield. For large values of W/b, the structure is more easily 
seen as a ridged waveguide loaded with dielectric. When 
W/b =1, a dielectric slab loaded waveguide exists and finally 
for W/b = 1 and = 1, the structure takes on the configura- 
tion of empty rectangular waveguide. 

The most commonly used transmission line to date in 
microwave integrated circuits is the microstrip line. This 
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Figure 1. 3-Dimensional View of Fin-line Structure 
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structure when used at high frequencies (millimeter wave 
range) results in problems which include radiation loss, 
spurious coupling, dispersion and higher mode propagation. 

In fin-line structures, metal fins are printed on a 
dielectric substrate which bridges the broad walls of a 
rectangular waveguide. In effect, the line is a printed 
ridged waveguide and can be designed to have a wider useful 
bandwidth than conventional line and can provide bandwidths 
in excess of an octave with less attenuation than microstrip. 
This adaption of ridged waveguide permits circuit elements 
to be fabricated by photo-etching at low cost and is compati- 
ble with thin-film hybrid techniques. 

In passive circuits, such as filters, the fins as illus- 
trated in Figure 1 may be directly grounded to the metal walls 
of the waveguide, and lumped elements, such as beam-lead 
capacitors, may be added. The gap between the fins can be 
varied along the longitudinal axis to provide low-cost circuit 
elements. When fin-lines are constructed as illustrated by 
Figure 2, the waveguide is parted along a plane where the 
current flow is parallel to the break, as in a common slotted 
line and the dielectric material extends through the separation 
in the broad wall of the shield. The upper fin is insulated 
from the housing at dc by a dielectric gasket. By making the 
thickness of the broad walls a quarter wavelength, a short 
circuit will appear between the fin and the inner wall of the 
shield at RF. Thus a practical fin-line will have the same 
electrical characteristic as the idealized structure of 
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Figure 2. Integrated Fin-line Semiconductor Mount 
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Figure 1. The lower fin is grounded directly by a metal 
gasket to provide a dc return when solid state devices are 
mounted on the fins. 

B . RELATED WORK 

Fin-line was first proposed in 1972 by Meier [Refs. 2 
and 3] and has since received considerable attention by a 
number of other investigators. Meier's work was primarily 
experimental in nature and the wavelength and impedance 
calculations were approximated. Previous analytical work on 
fin-line has been based on various methods [Refs. 8-11] . 

In this study, dispersion of different fin-line geometries 
was obtained by applying Galerkin's method to the electric 
fields in the slot interface [Ref. 8] . Another approach used 
the Transmission Line Matrix (TLM) technique yielding the 
dispersion characteristic of the fundamental and higher order 
modes of propagation [Ref. 9]. 

C. STATEMENT OF OBJECTIVES 

The Spectral Domain technique is used here in the analy- 
sis of the fin-line structure shown in Figure 1. This tech- 
nicue was first suggested by Itoh and Mittra [Ref. 12] and 
has since been applied to the analysis of other structures 
[Refs. 13-15]. 

One of the advantages of this approach is that it is 
numerically more efficient than the conventional methods that 
work directly in the space domain. This is due primarily to 
the fact that the process of Fourier transformation of the 



15 



« 



n 




coupled integral equations in the space domain yields a pair 
of algebraic equations in the transform domain that are 
relatively easier to handle. Another important advantage is 
that the Green's function takes a much simpler form in the 
transform domain, as compared to the space domain where no 
convenient form of the Green's function is known to exist. 

The Spectral Domain Transform method is first applied 
here in the formulation of the dispersion characteristic of 
the fin-line structure of Figure 1. During the application, 
a new matrix approach to the implementation of this method 
is presented. The characteristic impedance is next formu- 
lated using this technique in terms of the dispersion char- 
acteristic. Following the theroetical analysis, an explana- 
tion of the computer program used in determination of the 
wavelength and characteristic impedance is presented. Numeri- 
cal results are then compared with known data for ridged 
waveguide, slab loaded waveguide and slot line. These com- 
parisons establish the accuracy of the numerical results and 
illustrate the applicability of this method for the full 
range of structure parameters. Finally, several families 
of design curves are presented for a practical choice of 
f in- line parameters . 
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II. THEORETICAL ANALYSIS OF FIN-LINE 



A. FIELD AND BOUNDARY CONDITIONS 

Fin-line supports a hybrid field and it is known that all 

hybrid field components can be obtained from the superposition 

of TE and TM modes which are related to the two scalar poten- 

0 h 

tial functions 4i (x,y) and (p (x,y) , where the superscripts e 
and h denote electric and magnetic, respectively [Ref. 16]. 

The axial components of TM and TE are then 



T. 1 2 ^ e , . rz 

4 , (x,y) e 



( 1 ) 



TT 1 2 , h , . rz 

= k^ (j> (x,y) e 



( 2 ) 



2 2 2 2 . 
where k . = k. + r and k. =o))j.e., i=l, 2, 3 for each 

Cl 1 111 

of the three regions defined in Figure (1) . 

Through Maxwell's curl equations the transverse field 
components are then determined by these axial components and 
can be given as 



E = (r 

X 


dx 


- j 


34) , rz 
“ “ 3y > ® 


(3) 


E = (r 

y 


34)® 

9y 


+ j 


34) , rz 
“ ” 5x > ® 


(4) 



H 



X 






rz 



(5) 
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( 6 ) 



H 



(r 



8y 



- j 0) £ 



8 i|) 

3x 



rz 



since the lossless case is assumed, the propagation constant 

2 2 2 

is r = ±j6 and k . = k. - 6 . We will also assume that 

Cl 1 

= £, = e and e- = e_e„. 

1 3 o 2 o r 

Applying boundary conditions at the walls in region 1, 
tangential field components must be zero, so it follows that 
at y = hj^ + D 

E , (x,h- +D, z) = 0 

Z 1 1 

Exi (x,hi+D, z) = 0 

At the interface between region 
components must be continuous. 

At y = D 



3H ^ 

-^(x,h^+D,z) = 0 (7) 

9H , 

■ Yy -- (x,hj^+D, z) = 0. (8) 

1 and 2, tangential field 



E , (x,D,z) = E „(x,D,z) (9) 

z 1 z z 

E^l(c,D,z) = ^ (10) 

Also the electric fields at y = D will exist only in the 
slot and can be expressed as. 
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E^, (x,D, z) 
zl 



0 I X I ^ W/ 2 

e^(x)e^^ |x| < W/2 



( 11 ) 



E , (x,D,z) 
xl 



0 I^I 1 W/2 

e (x) jx] <W/2 



( 12 ) 



Similarly, tangential magnetic fields must be discontinuous 
by corresponding surface current densities. 



H , (x,D,z) - H ,(x,D,z) 
z 1 z z 



j (x) e^^ |xl ^ W/2 

X 

(13) 

0 |xl < W/2 



H^^(x,D,z) - 



j (x) e^^ I X I ^ W/2 
0 |xl < W/2 



(14) 



The tangential field components at the interface between 
region 2 and 3 must also be continuous. At y = 0 



E^, (x,0 , z) 
Z2 



E^2 (x, 0 ,z) 



(15) 



E^2 ( X , 0 , z ) 



Ex3 (x, 0,z) 



(16) 
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(x,0,z) 



H _ (x,0, z) 
Z J 



(17) 



Once again at the shield wall in region 3, the tangential^ 
field components must be zero, therefore at y = -h 2 

3H 

E^^(x,-h 2 ,z) = 0 -- g - y --(x,-h 2 f z) = 0 (19) 

8H ^ 

E^ 2 (^'~h 2 /Z) = 0 -- g^ --(x,-h 2 f z) = 0 . (20) 

The final boundary conditions occur at x = ±b/2 where the 
tangential components must be zero in all regions. At 
X = ±b/2 

E^^(±b/2,y,z) = 0 (21) 

Exi(±b/2,y,z) = 0 . (22) 

B. SPECTRAL DOMAIN APPROACH TO DISPERSION CHARACTERISTIC 

The approach used in the analysis of the dispersion 
characteristic is basically a modification of Galerkin's 
approach adapted for application in the Fourier transform 
or spectral domain. 

0 K 

The scalar potential functions (}) and <j) satisfy the 
Helmholtz equations in the three spatial regions , thus 
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I 



i 



0 



(23) 



't * 



Cl 



2 

where denotes the two-dimensional Laplacian operator in the 
transverse direction. These equations are second-order differ- 
ential equations for the unknown potential functions which 
may be Fourier transformed with respect to x to obtain ordinary 
differential equations via the transformation defined as 



~ j a X 

F{(f)(x,y)} = $(a^,y) = / <j> (x,y) e ^ dx . (24) 

— 00 



The Helmholtz equations of (23) are transformed into 



3 ^i(“j^/y) 
3y^ 



(a^ - k^. ) <}>. (a ,y) 

n Cl 1 n ^ 



(25) 



0 

for the electric field, <J> , and the magnetic field, <f> . 

The solutions to these ordinary differential equations 
for the three regions after applying the boundary conditions 
at y = h^ + D and y = -h 2 (equations (7), (8), (19) and (20)) 

can be written as 



'J>^(ctnfy) = A®(a^) sinh Y 2 ^(D+h^-y) 



(26a) 



^^(a^,y) = 



B®(a^) sinh Y 2 y+C®(a^) cosh Y 2 y (26b) 
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(26c) 



(a 

3 n 


^y) 


II 

D 

(D 

P 


sinh 


Y3(h2+Y) 


*l<“n 


/Y) 




cosh 


Yi(D+hi-y) 


(a 

2 n 


/Y) 


= B^a^) 


sinh 




*3<“n 


/Y) 


Q 

II 


cosh 


Y3(h2+Y) 



(26d) 



(26e) 



(26f) 



where 



2 , 2 
a - k . 
n Cl 



2 2 2 
a + B - kf 
n 1 



(27) 



and 



n2iT/b 



even 



(28a) 



a 

n 



(2n-l) TT/b 



odd 



(28b) 



2 

It is important to observe at this point that Yj_ inay be 

less than zero in any of the three regions of the structure 

under certain conditions. Therefore three sets of solutions 

can exist for equations (26a) to (26f) . When = 0 and k^ 

approaches k^ (where k^ is the wave number for free space) 

2 

B is less than k. and so y- < 0. Under this condition the 

1 1 

hyperbolic functions in all three regions are replaced by 

2 2 

trigonometric functions. If k < B < k^, then y, and y_ 

0 2 13 
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2 

are greater than zero and Y 2 ^ some values of and 

the trigonometric functions replace the hyperbolic functions 

in the spatial region 2 only. This suggests that the nature 

of the field is dependent upon the values of the transform 

2 

variable, a . The last solution occurs when y. >0 for which 
n 1 

2 

the solutions are given above. For conditions when y^ < 0, 

II **2 2 

replaces Yj such that (Yj_) = and cosh Yj^y is replaced 

II II 

by cos Yj^y and sinh Yj^y is replaced by j sin Yj^y* 

The next step in obtaining a solution to the differential 

equations is the application of the continuity or boundary 

conditions at the interfaces between the three spatial regions, 

namely at y = 0 and y = D. These conditions are given in 

equations (9) through (18) . Through their application, the 

0 h 

eight coefficients A through D may be related to each 
other, to the fin surface current j (x) and the slot field 
e(x). The resulting set of linear equations (really a set of 
transformed boundary equations) may be written in a matrix form 
as follows: 










0 


e 






B 




0 


• 




• 


• 




• 


• 




• 






Ex 
, Ez 



(29a) 
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I 



t 



I 

i 






(29b) 



A® 




0 






0 


• 




. > 


• 


II 


1 • 


• 




• 


. • 




J 

X 






J 

z 



where E. (a ) is the transform of the slot field and J. (a ) 
i n in 

is the transform of the surface current. The matrices [M^] 

and [Mj] are square 8x8 matrices and differ only in the 

last two rows. Details of the application of the continuity 

conditions and the formation of the [M^] and [Mj] matrices 

are shown in Appendix A, with a listing of the matrix elements 

2 2 

under the conditions y. > 0 and y- *^0 ^nd normalized with 
respect to D. 

To solve the m.atrix equations (29) , the unknown coeffi- 
• 0 h 

cients A through D are determined by 



A®' 

• 

• 






0 

0 








E 

X 

E 

z 



where [M^] ^ is the inverse matrix of [M^] . Now substituting 
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equation (30) into (29b) results in the single matrix 
equation 





0 




0 




0 




0 


[Mj] [M^] ^ 


• 

• 


= 


• 

• 




Ex 




Jx 




Ez 




. Jz 



Using the four lower right hand corner elements of the 
matrix [Mj] [M^] (31) is converted to 



■ G, (a , B ) (a , B 

In Z n 

G -5 (a / 6 ) G . (a / 3 

3 n 4 ' n 





E (a ) 




J 


(a ) 




X n 




X 


n 








. J 


(ci^) 



(32) 



where the elements of the matrix [G] are the Fourier transforms 
of the components of the dyadic Green's function for this 
structure . 

A solution to equation (32) is obtained by using the 
Method of Moments [Ref. 17] . The inner product is defined as 



oo 

<f(a^), g(a^)> = I f(a^) g*(a^^) . (33) 

n=-oo 
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Galerkin's method is used here and thus the Fourier trans- 
forms for the electric field components, E and E will be 

^ z 

used for the weighting functions, W(a^). Taking the inner 

product of equations (32) with 

W-(a ) = E (a ) yields 
2 n z n 



Following from Parseval's theorem, the right-hand side of 
these equations are zero because of the orthogonality of 
the electric field and surface current at y = D. 

Up to this point, the formulation of the problem is exact 
since no approximations have been made. The electric field 
components and their transforms although unknown can be 
expanded in a set of basis functions. An investigation of 
various one-term approximations made by Knorr and Kuchler 
[Ref. 18] suggests good accuracy can be achieved if the 
electric field distribution between the fins is approximated 
by 



<S^x'^x>-^<S^z'^x> = 0 



(34) 



<G-E ,E >+<G.E ,E.,> = 0 . 

^ 3 X z-^ 4 X z ^ 



1 



|x| <_ W/2 




0 



elsewhere 



( 35 ) 




0 
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for the dominant mode. This choice (shown in Figure 3 for 
A = 1) has been shown to give accurate results for slotlines 
with W/D ^ 2 [Ref. 18] and it is exact for dielectric slab- 
loaded waveguide (W/b = 1). Assuming that e^(x) = 0 reduces 
equation (34) to 

oo 

n=-“ 

where E (a ) is the transform of the electric field distri- 
X n 

bution 



E (a ) 

X n 



W/ 2 j a X 

/ e (x) e ^ dx 

-W/2 



AW 



sin (a W/2) 
n 

w/2 . 

n 

(37) 



So, one can see that (36) is dependent on the ratio of X/X’ 

or effective dielectric constant, £ , defined by 

^eff 



= (X/X')^ = 6^/k^ . 

^eff 



The dispersion characteristic of the dominant fin-line 
mode can now be found by varying 6, the propagation constant, 
such that equation (36) is satisfied for a given set of 
physical parameters at a desired frequency of operation. 
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Figure 3. Assumed Electric Field Component in Slot 
in x-direction Versus x for Fin-line 
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C. CHARACTERISTIC IMPEDANCE IN TERMS OF DISPERSION 

CHARACTERISTICS 

The definition of the characteristic impedance for an 
ideal TEM transmission line is uniquely given by static 
quantities. Since fin-line is not a TEM transmission line 
nor does it support a pure TE or TM waveguide mode but rather 
supports hybrid modes, no unique definition of the character- 
istic impedance can be found. 

One possible choice is to define it as 



z 

o 



V 



2 P 



avg 



(38) 



where is the slot voltage defined as 



V = - [e - di . 

° slot 

Assuming the electric field shown in Figure 3 and (35), 



V 

o 



W/2 

/ A dx 

-W/2 



1 



where A was arbitrarily selected as 1/W so that AW = 1. 

Pavg ^he time-averaged power flow in the fin-line struc- 

ture which is given by 
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( 39 ) 



P = ^ Re / / E X H*- a da 

avg 2 z 

= Re / / (E H * - E H dx dy . 

2 X y y X 

Substituting equations (3), (4), (5) and (6) into (39) 

results in 



avg 



I 



b/2 



D+h, 



-b/2 



-h. 






84 ) 



[(jef^- jcou^) (-je^+ 



ay 



8y 



3x 



- (jSff-+jo.p 



84)^, , . „84>^ . 84)^ 



•) ] dxdy , 



(40) 



This formulation applies to the time-averaged power flow 

in the spatial domain. Through the use of Parseval ' s theorem, 

it is possible to transform this expression into the spectral 

domain where the scalar potential functions are known. The 

integration with respect to x is transformed via Parseval 's 

relation into an infinite summation. Applying this theorem 
0 h 

where 4> (x,y) and 4> (x,y) are real yields 



b/2 



D+h. 



-b/2 



-h. 



4>® (x,y) 4)^ (x,y) dxdy 



= K E / 



D+h, 



-h. 



4'^ (oij^fY) 
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Therefore, equation (40) is transformed into 



avg 2b 



Re I 



n=-“ -h. 



^•"^1 e 

r r . 2|^ei2 „ 2. h,2 „ , 3<I> i2 

j {-6o)ea^|<I> I -Bdjvcij^ I $ | -6we|g-^| 



„ |8$^|2 ,„2 r^e a($^)* . a<l>® /*h,*T 

- - 36 Oj,[* — 3j— +3^(*)*1 



^ .,2 r*h 9('t'®)* ^ , 

+ j k a [ <I> + -r ( <I> ) * ] } dv . 

dy dy 



n 



(41) 



This expression must be evaluated in each of the three 
regions of the fin-line shown in Figure 1. Therefore the 
power flow may be expressed by 



P = Pt + P_ + P_ 

avg 123 



where P^^, ^3 defined as 



D+h, 



?1 = ^ I / {f (<&®,<&^) } dy 



n=-oo D 



(42a) 



Pj = ^ Re J J (f(*|,42))dy 



n=-<» 0 



(42b) 



‘ I 



0 , 

/ {f ($ 2 /^ 3 ) ^dy 



n=-oo - 






(42c) 
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The simple functional dependence on the variable y allows 

one to perform the integration in equations (42a) -(42c) 

analytically (see Appendix B) . This leaves P as the 

avg 

summation of the power flow in each region. 



avg 



2b 



n=-co 



(P^ + + Pj) 



(43) 



Now that can be determined by the application of the 

dispersion characteristic or propagation constant determined 
earlier and that the slot voltage is known, the character- 
istic impedance can be computed. The details of these lengthy 
but straightforward algebraic manipulations are shown in 
Appendix B . 
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III. COMPUTER PROGRAMMING 



In Section II the theoretical background of the analysis 
was provided. This section concerns itself with the computer 
programming technique used for fin-line analysis. The program 
was designed for the determination of the wavelength and the 
characteristic impedance of a particular fin-line structure 
over a range of frequencies but due to the nature of the 
structure various other configurations may be considered by 
the proper selection of input parameters. Dielectric slab 
loaded and ridged waveguide along with rectangular waveguide 
may also be analyzed. Before the program is described some 
general considerations which lead to its development are 
necessary for understanding of the numerical methods used. 

A. REMARKS ON NUMERICAL ANALYSIS 

The computation of the characteristic impedance is based 
upon the solution to the dispersion characteristic problem 
for the transmission line under consideration. In other 
words, the wavelength ratio, X'/\, or the wave propagation 
constant, 6, must be known before any other investigations can 
be started since only in this case are the scalar potential 
functions in the transform domain known. Hence, in the analy- 
sis of the various structures, the starting point is a solu- 
tion for 6/ the propagation constant of equation (36). 
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In general a 3 is sought such that 



I G, (a^,B) |E(a„) 1^ = 0 (44) 

^ in ' n ' 

n=-oo 

where Gj^(a^,S) is a component of the transformed dyadic 
Green's function and is the approximation to the 

Fourier transform of the electric field in the slot in the 
plane y = D. A first order approximation is assiimed here 
which eliminated the other terms of equation (34) . 

It is obvious that an analytic solution for 3 cannot be 
obtained due to the algebraic complexity leading to the quan- 
tity Gj^(a^,3) as shown in Appendix A. Thus a numerical 
solution is mandatory and must be obtained by computer 
programming . 

It is shown in Section II that the transforms of the dyadic 
Green's function components were obtained by mathematical 
manipulation of matrices. The matrix elements are either 
real or imaginary so that in general a complex computation on 
the computer is necessary. Furthermore, it is shown that the 
elements change from hyperbolic to trigonometric functions. 

If 

For examaple sinh ^ 2 ^ becomes jsin ¥ 2 '^* it found that 

these quantities are not only pure real or pure imaginary 

but change from real to imaginary or vice versa for the 
2 2 

cases Y • >0 and y • < 0 . 

1 1 

Previously, equation (32) was obtained by lengthy alge- 
braic manipulation of the boundary equations. In the approach 
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described here the matrices (29a) and (29b) are defined 
directly from the boundary equations and equation (32) is 
arrived at by numerical computation. There is some sacri- 
fice in numerical efficiency in this approach but the formu- 
lation of the problem is straightforward. 

Since in equation (44) the absolute value squared of the 
assumed electric field is taken, this part of the series will 
be in general a real and even function of a^. Preliminary 
numerical investigations by Kuchler [Ref. 19] of the coeffi- 
cients of these series indicated an even distribution with 
respect to the variable a^, also exists. Therefore, all terms 
of the series (44) are even functions which allows the saving 
of computation time by summing over only half the interval. 

A desired solution for equation (44) can now be found by 
assigning an arbitrary numerical value for the wavelength 
ratio and then iteratively changing this until the series con- 
verges to zero or more precisely a small enough value so 
that a prescribed accuracy for the wavelength ratio is realized. 

The next task is the determination of the time average 
power flow. Here again complex arithmetic expressions are 
involved as seen in Appendix B. 

Dielectric substrates are commercially available for 
various thicknesses D so that in general the line parameters 
should be computed for each thickness and frequency. A more 
generalized method of analysis is simply to scale the various 
geometric quantities involved with respect to this thickness 
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to .obtain more versatile parameters for a varying normalized 
frequency D/X . The resulting design curves may then be 
applied to the structures scaled to different operating 
frequencies . 

B. COMPUTER PROGRAM DESCRIPTION FOR FIN-LINE CONFIGURED 

TRANSMISSION LINE 

A computer program which provides a numerical analysis of 
fin-line is shown in Appendix C with the main program called 
'FIN-LINE' and its accompanying subroutines 'CALC' and "CMTRIN, ' 
all written in the FORTRAN IV programming language. 

The program 'FIN-LINE' starts by reading the input data 
which includes the structure dimensions, the dielectric con- 
stant, the initial values of the slot width and frequency 
and the iteration steps. After the input data is read, 
initialization of computation loops occur for varying slot 
widths and normalized frequencies. Next, the program determines 
the type of structure that exists and whether or not the con- 
figuration is operating in the dominant mode. Operation at 
higher order empty waveguide modes are allowed but program 
termination occurs at the onset of surface waves. 

The program proceeds with preliminary computations used 
in the determination of the dispersion characteristic. The 
first major computation of the program that takes place is 
the finding of the wavelength ratio X '/X which satisfies 
equation (36) . The necessary coefficients for determining 
the series are provided by the subroutine 'CALC. An initial 
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e is the lower 
r 



value of \'/\ = l//e^ is selected where 

limit of X'/A/ that is when A' = X., this wavelength of the 

a 

dielectric material. The results from the initial computation 

are used to select a new value for X'/X through a bisectional 

zero finding method to improve the results and the procedure 

starts over again. The computation is terminated when the 

correct value for X'/X is determined which satisfies equation 

(36) or when the absolute value of the difference between the 

-4 

old and the new value for X'/X is less than 5*10 and thus 
produces a three digit accuracy after the decimal point. 

Using the determined value of the wavelength ratio, the 
program continues with the computation of the characteristic 
impedance. The coefficients for this computation (equation 40) 
are provided again by the subroutine 'CALC'. 

After completion of the dispersion and impedance compu- 
tations, the program prints these values and restarts with 
new values for the normalized frequency and/or slot width. 

Figure 4 presents a flow diagram for the main program 
'FIN-LINE'. All the variable names have been chosen such 
that the statements can be readily compared with the respec- 
tive equations of the appropriate appendices. 

The subroutine 'CALC provides the coefficients and matrix 
elements as well as an extensive amount of the computation 
used in the determination of both the dispersion and impedance. 
The entry variable is X/X' for both characteristics. A flow 
diagram for the function 'CALC is shown in Figure 5. 
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Figure 4a. Flow Diagram for Main Program 'FIN-LINE' 
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Figure 4b. Flow Diagram for Main Program 'FI>7-LIME' 
(Continued) 
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Figure 4c. Flow Diagram for Main Program 'FIN-LINE' 
(Continued) 
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Figure 4d. Flow Diagram for Main Program 'FIN-LINE' 
(Continued) 
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Figure 4e. Flow Diagram for Main Program 'FIN-LINE' 
(Continued 
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Figure 4f. Flow Diagram for Main Program 'FIN-LINE' 
(Continued) 
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Figure 5a. Flow Diagram of Subroutine 'CALC' 
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Figure 5b. Flow Diagram of Subroutine 'CALC' 
(Continued) 
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Figure 5c. Flow Diagram of Subroutine 'CALC 
(Continued) 
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Figure 5d. Flow Diagram of a Subroutine 'CALC' 
(Continued) 
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Figure 5e. Flow Diagram of Subroutine ’CALC’ 
(Continued) 
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Figure 5f. Flow Diagram of Subroutine 'CALC' 
(Continued) 
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Figure 5g. Flow Diagram of Subroutine 'CALC 
(Continued) 
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Figure 5h. Flow Diagram of Subroutine 'CALC' 
(Continued) 
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Figure 5i. Flow Diagram of Subroutine 'CALC 
(Continued) 
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Figure 5j . Flow Diagram of Subroutine 'CALC 
(Continued) 
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Figure 5k. Flow Diagram of Subroutine 'CALC' 
(Continued) 
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Once the subroutine is entered from the main program, 

variables used throughout are computed. Following the initial 

computations, a decision is made as to whether hyperbolic or 

trigonometric functions are to be used in the formation of 

the matrix elements shown in Appendix A. This decision, of 

2 

course, is made by determining whether is less than or 

2 

greater than zero. If > 0/ hyperbolic functions apply 
2 

and if < 0/ trigonometric replace the hyperbolic functions. 
After the decision making and the formation of the [ME] and 
[MJ] matrices, the [ME] matrix is inverted and the matrix 
elements of concern in the [MJ] and [ME] ^ matrices are multi- 
plied together for the formation of element Gll of the dyadic 
Green's functions. 

The transform, of the assumed slot electric field distribu- 
tion is determined. At this point the decision, as to whether 
the dispersion or the impedance is being determined, is made. 
If the dispersion computation is in process, the subroutine 
computes the terms of equation (36) and returns to the main 
program for a new selection of A'/X. 

If the characteristic impedance computation is being per- 
formed, the program continues with the determination of the 

• . s h 

coefficient matrix elements A through D . Once again the 

decision for the use of trigonometric or hyperbolic function 

is made followed by additional variable computations needed 

in determination of power. 
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The power flow in each region is determined for 
greater or less than zero, depending on the conditions, and 
totaled. This procedure is followed until the difference 
between two consecutive values of the impedance is less than 
0.1% at which time return to the main program occurs. 

The subroutine 'CMTRIN' takes the complex [ME] matrix 
and inverts it into the [ME] ^ matrix destroying [ME] in the 
process. The Gaussian Elimination Method with column pivoting 
is used in determination of the inverse of the matrix. Since 
this subroutine is small and performs only the one function, 
it has not been flow charted but is listed in Appendix C. 

C. PROGRAM LIMITATIONS 

Infinite summations such as those which appear in equa- 
tions (36) and (41), can not be performed on a digital com- 
puter which deals only with finite numbers. Therefore infinity 
must be replaced by a finite approximation. In the case of 
equation (36), an analysis of ^ versus n was con- 

ducted with results shown in Figure 6. As shown in this figure, 
when n reaches a value of about 10 the expression is down by 
at least a factor of 10 ^ relative to its largest value close 
to the origin. The value n = 20 is sufficiently large to 
be considered infinite for all practical purposes. Therefore 
equation (36) becomes 



00 

I 1 E 

n=0 ^ 



Gj^(O) 



20 

I 

n=l 



G, (a , 6) 
1 n 



E (a ) 
X n 



(45) 
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Figure 6. | ^ vs n with = 2.2 and 

D/A = 0.02 for Five Arbitrary Values 
of X'/x 
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After establishing this limit on the summation, computer 
errors (overflows) occurred at various points in the summa- 
tion. An investigation revealed that the problem originated 
in the dyadic Green's function component when became 

too large. Varying the size of h^ produced overflows only 
when h was increased beyond a certain point; a situation 
where the walls around the central slotline-like structure 
became far removed from the field. Limits were able to be 
placed on such that if its value became greater than some 

value X in this case X = 80, then Yj_hj was set equal to X 
because of the limited interaction between the walls and the 
field. When the field-wall interaction becomes a factor, 
the overflow problem no longer exists. Limitations of this 
kind are found in the computations of the wavelength as well 
as the impedance. 

The high valued arguments of the hyperbolic functions just 

discussed cause additional problems when determining the 
0 h 

coefficients A through D in the power flow computations. 

In the inversion of the matrix [ME] , these large values become 

small ones and when squared produce underflows in the computer 

system since the magnitude of the elements were less than 

10 Limitations were placed only on the coefficients A^, 

0 

D and D which caused the underflow problem. The coefficients 

— 38 

were limited to magnitudes of 5*10 . This problem was 

compounded by the fact that the coefficients are complex and 
they may also take on either plus or minus signs. 
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During the investigation of the overflow problem, it 
was noted that when W/b = 1, E =0 for n > 0 and so G, (0) 
was the only important term in the computation. Therefore 
whenever W/b = 1 only the n = 0 term is computed in both 
the wavelength and impedance computations. 

In the characteristic impedance equation (41) , the 

coefficients in this infinite series were found to decay 

rapidly so that a finite approximation yields good results. 

Figure 7 shows the characteristic impedance as a function of 

N, the number of terms in the truncated series, for various 

slot widths. Termination of the impedance computation occurs 

when the difference between two consecutive values of in 

o 

the iteration loop falls below 0.1% or when the upper limit 
of N = 50 is reached. A minimum of seven iterations is set 
for the larger values of W/b except for the situation when 
W/b =1. In this case, only the N = 0 term is non-zero. 

D. INPUT/OUTPUT DATA 

Proper input data formatting is essential for any success- 
ful computer program. Input data, similar to that shown in 
Tables I and II are required with the program 'FIN-LINE'. A 
minimum of three data cards is necessary to make up the 
input data deck. The first card specifies the dimensions of 
the shielding and the permittivity of the dielectric material 
(b/D, hj^/D, h 2 /D and e^) . This information is entered in a 
floating point format with the column requirements specified 
in the program listing in Appendix C. Note that for a ridged 
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Figure 7. Characteristic Impedance, Z^, vs. N, the 
number of terms in the truncated series, 
with e =2.2 and D/X = 0.0212 for a WR(19) 
shield^ 
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Table I. Sample Input Data for a Ridged Waveguide 
Structure with WR(19) Shield 



18.8 

3 


18.8 


17.8 


2.2 




0.02 


0. 0170 


15 


0. 0 


0. 0021 


0.5 
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15 


0.0 


0. 04 


1.0 


0.01 


15 


0.0 


0. 01 
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Table II. Sample Input Data for Fin-line Structure 
with WR(19) Shield 
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waveguide, the dielectric substrate thickness is zero. This 
would result in infinite values for the first three values 
of the input data (b/D, h^^/D and ^ 2 / 0 ) . To eliminate this 
problem select any thickness greater than zero, but let its 
dielectric constant equal one. 

The second card specifies the number of initial value 
data cards that follow with nine being the maximum allowed. 
Sample input of Table I indicates that two cards will follow 
the second card and sample input of Table II shows three data 
cards follow. 

The third card needed for the input contains initial 
values for the computations that follow. Six pieces of informa- 
tion are supplied by this card. The first two pieces of 
information contain the initial values of W/b and D/A. Next 
a two digit number specifies the number of different values 
of W/b (first digit) and D/A (second digit) that will be com- 
puted. On the last card of Table I, the n\imber 24 indicates 
that the wavelength and impedance will be determined for two 
different values of W/b and for four different values of D/A 
for each W/b. The last two values are the incremental values 
of W/b and D/A, respectively. The value of 0.5 will be used 
to increment the initial value of W/b so that the second value 
will be W/b = 1. The value 0.0021 will be used to increment 
the initial value of D/A so that the four values of D/A are 
0.0212, 0.0233, 0.0254 and 0.0275. 

The data shown in Tables III and IV is the computer output 
from execution of the program listed in Appendix C with the 
input data of Table I and Table II, respectively. 
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SLOT -TRANSMISSION LIN? ON A DIELECTRIC SUBSTRATE 
WITH SHIELDED WALLS 
(FIN-LINE) 



DIMENSIONS: B/D=18.80 Hl/D=18.80 H2/0=17.80 5PSR= 1.0 



RIDGED WAVEGUIDE STRUCTURE EXISTS 
FOR THE FOLLOWING VALUES OF W/B 



W/B= 0.100 


D/L = 


0.0170 


LP/L= 1.1559 


Z0= 215.7 


OHMS 


W/B= 0.100 


D/L = 


0.0191 


LP/L= 1.1176 


Z0= 208.5 


OHMS 


W/B= 0.100 


C/L = 


0.0212 


LP/L= 1.0918 


Z0= 203.8 


OHMS 


W/B= C.lOO 


D/L = 


0.0233 


LP/L= 1.0738 


Z0= 200.5 


OHMS 


W/B= 0.100 


0/L = 


0.0254 


LP/L= 1.0613 


Z0= 198.1 


OHMS 


W/B= 0.500 


C/L = 


0.0212 


LP/L= 1.1988 


Z0= 382.6 


OHMS 


W/B= 0.500 


D/L = 


0.0233 


LP/L= 1.1559 


Z0= 369.1 


OHMS 


W/B= 0.500 


C/L = 


0.0254 


LP/L= 1.1262 


Z0= 359.6 


OHMS 


TE 


-20 EMPTY GUIDE MODE PROPAGATES 
VALUES OF D/L GREATER THAN 0.027 


FOR 




W/3= 0.500 


0/L= 


0.0275 


LP/L= 1.1051 


Z0= 352.5 


OHMS 




RECTANGULAR WAVEGUIDE STRUCTURE 
EXISTS FOR W/B= 1.0 






W/B= 1.000 


D/L = 


0.C212 


LP/L= 1.2840 


Z0= 484.1 


OHMS 


W/B= 1.000 


D/L = 


0.0233 


LP/L= 1.2176 


Z0= 459.2 


OHMS 


W/B= 1.000 


0/L = 


0.0254 


LP/L= 1.1738 


Z0= 442.4 


OHMS 


W/B= 1.000 


D/L = 


0.0275 


LP/L= 1.1426 


Z0= 430.6 


OHMS 



Table III. Output Data Obtained Using Sample 
Input Data from Table I 
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SLOT -TRANSMISSION LIN5 ON A 0I5LECTR I C ' SUBSTRATE 
WITH SHIELDED WALLS 
(FIN-LINE) 



DIMENSIONS: B/D=18.80 Hl/D=18.80 H2/0=17.80 EPSR= 2.2 



W/B= 0.020 


D/L= 0.0170 


LP/L= 0.9183 


Z9 = 


138.6 


OHMS 


W/3= 0.020 


D/L= 0.0191 


LP/L= 0.9035 


Z0 = 


138.1 


OHMS 


W/B= 0.020 


D/L= 0.0212 


LP/L= 0.8933 


zo= 


138.2 


OHMS 


W/B= 0.020 


D/L= 0.0233 


LP/L= 0.8855 


zo= 


138.9 


OHMS 


W/B= 0.020 


D/L= 0.0254 


LP/L= 0.8793 


zo= 


139 .8 


OHMS 


W/B= 0.500 


C/L= 0.0170 


LP/L= 1.2457 


zo= 


421.2 


OHMS 


TE- 


-20 EMPTY GUIDE 
VALUES OF D/L GR 


MCDE PROPAGATES FOR 
EATER THAN 0.027 






W/B= 0.500 


D/L= 0.0570 


LP/L= 0.9465 


zo= 


501.2 


OHMS 


W/B= 0.500 


D/L= 0.0970 


LP/L= 0.9043 


zo= 


700 .0 


OHMS 


W/B= 0.500 


0/L= 0.1370 


LP/L= 0.8683 


zo= 


776.9 


OHMS 


SURFACE WAVES PROPAGATE FOR VALUES OF 
GREATER THi^N 0.169; RUN TERMINATED 


C/L 






DIELECTRIC SLAB-LOADED WAVEGUIDE 
STRUCTURE EXISTS FOR W/8= 1.0 






CUTOFF FREQUENCY OF THE DOMINANT 
OF EMPTY GUIDE IS C/L= 0.013 


MODE 






W/B= 1.000 


D/L= 0.0100 


Lp/L= 


Z0 = 






W/B= 1.000 


D/L= 0.0200 


LP/L= 1.2660 


Z0 = 


508.3 


OHMS 


W/B= 1.000 


0/L= 0.0300 


LP/L= 1.0707 


Z0= 


464 .9 


OHMS 


W/B= 1.000 


D/L= 0.0400 


LP/L= 1.0191 


Z0 = 


493.8 


OHMS 


W/B= 1.000 


C/L= 0.0500 


LP/L= 0.9957 


zo= 


553.6 


OHMS 



Table IV. Output Data Obtained Using Sample 
Sample Input Data from Table II 



64 



In addition to the tabulated results of the computations 
for the wavelength and impedance, message statements are also 
printed under certain conditions. All the messages that can 
be generated by the program are found in the output of Tables 
III and IV. 

Since various structures may take form by varying the 
input parameters, description of the type of structure that 
exists will be printed. For ridged waveguide, the permittivity, 
must be one and W/b < 1 whereas the structure becomes 
empty rectangular waveguide when W/b = 1 (shown in Table III) . 
With > 1 and W/b < 1, a fin-line structure is formed. Since 
this analysis is of fin-line structure, no output statement 
is printed. When W/b = 1 and > If a dielectric slab loaded 
waveguide structure exists (Table IV) . 

The analysis of the fin-line structure presented here 
concerns itself mainly with operation in the dominant TE^q 
limit mode. At the lower end of the dominant mode frequency 
range when D/X < D/X^, the cutoff frequency of the empty wave- 
guide, an output statement with this information is printed. 

The cutoff wavelength of the empty waveguide, X^, is equal to 
2a where a = hj^+h 2 +D. Since the structure takes other forms 

than just an empty waveguide, operation can exist below D/X . 

o 

with the presence of a dielectric slab in the fin-line and 

slab loaded waveguide operation exists below D/X . Therefore 

c 

analysis below the cutoff frequency of empty waveguide is 
possible but a lower limit is set. When the value of X'/X 
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exceeds 6, indicating the approach of cutoff, computations 
terminate without printing values of X'/X and for the 
values of D/X for which X'/X > 6. 

At the high end of the dominant mode frequency range when 
X = a, the TE 2 Q of the empty guide can propagate. Therefore 
a statement is printed when D/x ^ D/a. For a value of 
D/X ^ l/(4/7^), surface waves propagate on the metal-backed 
dielectric slab and normal operation is no longer possible 
so the run is terminated for that value of D/X. 

The total program length including all subroutines is 
48 K bytes of storage. The sample computer runs presented 
here with outputs shown in Tables III and IV took 4 and 5 
minutes, respectively which includes compilation time of 
about 20 to 30 seconds. 
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IV. NUMERICAL RESULTS AND COMPARISONS 



To check the accuracy of the numerical results generated 
by the computer program, comparisons were made with data 
available in the literature for various fin-line sub-structures. 
These structures vary from ridged waveguide at one end to 
dielectric slab loaded waveguide at the other. 

A. RIDGED WAVEGUIDE 

As mentioned earlier, the fin-line structure shown in 
Figure 1 becomes ridged waveguide with zero thickness ridges 
when = 1 and W/b < 1 or when the dielectric substrate 
thickness D is reduced to zero. 

The design of ridged waveguide was first treated by Cohn 
[Ref. 20] and later by Hopfer [Ref. 21] and Lagerlof [Ref. 

22] . The wavelength and impedance of the ridged waveguide 
structure can be given by the following expressions 
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The cutoff wavelength X^ was obtained from the curves pre- 
sented in [Ref. 21] for several values of W/b. Using this 
information and equation (46a), the wavelength ratio X'/^ was 
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determined. These calculated values were then compared with 
the results using the spectral domain method of the computer 
program discussed in Section III. Figure 8 shows the two 
results with agreement to within 0.5% where the discrete 
points are from equation (46a) . 

A similar comparison of the impedance was performed using 
equation (46b) where was obtained from curves found in 

[Ref. 22] and the numerical results of the program 'FIN-LINE I 
The results appear in Figure 9 where again very good agreement 
exists with a difference of only 2%. 

B. DIELECTRIC SLAB-LOADED WAVEGUIDE 

A variation of the fin-line structure of Figure 1 where 
W/b = 1 and > 1 results in a dielectric slab loaded 
rectangular waveguide. This structure has been studied in 
considerable detail by Vartanian, et al., [Ref. 23]. In their 
work, the dielectric slab was centered in the waveguide and 
they presented an analytical expression for the (voltage) 
impedance, Z (using the notation of [23]) at the center of 

p V 

the slab. The impedance computed in this analysis was speci- 
fied at the edge of the dielectric slab. Therefore a relation- 
ship between the two impedances must be defined before any 
comparison can be made. The impedance at the edge of the slab 
can be easily defined by the expression 
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Figure 8. Wavelength ratio, \'/^r vs. Frequency 
for a Ridged Waveguide. Ridges are 
Centered and Have Zero Thickness 
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Figure 9. Voltage Impedance Z vs. Frequency for 
a Ridged Waveguide.*^ Ridges are 
Centered and Have Zero Thickness 
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where is the field at the edge of the slab and 

is the field at the center of the slab. Thus, 

cos^(3£) (47b) 

o pv 2s 

where the various quantities are defined in [22] as 

= (2tt)^ [£ - (X/X')^] 

X 

s = c/2 

and C = D is the dielectric slab thickness. The characteristic 
impedance of a slab loaded waveguide was computed using 
equation (47b) (where Z^^ was obtained from Vartanian's 
curves) and the spectral domain method. The results of the 
computations are compared in Figure 10 where the discrete 
points are calculated from equation (47b) . It can be seen 
that the results agree very closely. The difference is less 
than 1%. 

The wavelength of the slab loaded waveguide may be 
determined analytically by the transverse resonance proce- 
dure [Ref. 16] which is found by writing the field expressions 
for an exact solution of the problem, applying boundary condi- 
tions at the air-dielectric interfaces and thereby obtaining 
the determinantal equation for the structure. The wavelength 
ratio is found by determining the value of X'/X which satisfies 
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Figure 10. Voltage Impedance Z vs. Frequency 
for a Slab Loaded wS(19) Waveguide. 
Slab is Centered 
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the determinantal equation. Figure 11 shows the variation 
of the wavelength ratio with dielectric thickness, D, for 
several values of obtained by the spectral domain method. 

The edge of the slab was placed at the mid-point of the broad 
wall of a WR(19) waveguide operating at 40 GHz. The dis- 
crete points are those values determined by the transverse 
resonance approach. Excellent agreement was obtained between 
the results of the two methods. The greatest difference between 
the different approaches was only 0.6%. 

C . SLOTLINE 

If W/D 2 and is sufficiently high for the fin-line 
structure of Figure 1, the presence of the shield will have 
little effect if the walls are sufficiently far removed from 
the slot. In this case the structure will behave like a slot- 
line. This behavior of fin-line is illustrated in Figure 12 
where the wavelength and characteristic impedance is plotted 
for = 20 and a slot width normalized to the dielectric 
thickness, W/D = 1. The discrete points represent the wave- 
length ( □ ) and the impedance ( 0 ) of a slotline obtained 
from [Ref. 24] for the same parameters and W/D as the fin- 
line. Agreement between the fin- line and the slotline is 
within 1% for the wavelength and 5% for the impedance. The 
impedance curve for a ridged waveguide structure with the 
same fin separation is plotted for reference. 
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Figure 11. Wavelength ratio, X'/X vs. Dielectric 
Thiclcness, for a Slab Loaded WR(19) 
Guide. Edge'^'of Slab is Centered 
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Figure 12. Wavelength Ratio, A '/A, and Characteristic 
Impedance, vs. Normalized Frequency D/A 
for Fin-line°and Slotline. Fin-line Has 
WR(19) Shield and Fins Centered 
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D . MEASUREMENTS 



Very little experimental data is available in the 
literature for any significant comparison. However , fin-line 
wavelength measureirients have been presented by Meier [Ref. 1] . 
The guide wavelength in the millimeter fin-line structure was 
measured across the band of 26.5-40 GHz by a sliding short- 
circuit technique. The fin-line used a D = 10 mil substrate 
with = 2.2. A WR(28) shield (a = .280", b = .140") was 
used along with a fin separation W/b = .128 (W = .018") . 

The results of this experimental work are plotted in Figure 
13. A computer run was made for a structure with the 
parameters reported by Meier and the results were compared 
with Meier's measured values. This comparison is illustrated 
in Figure 13 where an agreement to within 1% can be seen. 
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Figure 13. Computed and Measured Values of Wavelength 
Ratio X'/X vs. Frequency for Fin-line 
with WR(28) Shield. Fins are Centered 
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V. FIN-LINE DESIGN CURVES 



A problem exists with the preparation of design curves 
for structures like fin-line because of the mamber of inde- 
pendently variable parameters which describe the structure. 

This difficulty can be alleviated to some degree by practi- 
cal considerations. First, fin-lines are constructed with 
a shield that is compatible with the dimensions of standard 
rectangular waveguide for the millimeter wavebands. Above 
22 GHz, all the standard waveguide structures have aspect 
ratios, b/a = 0.5. Furthermore, the fins are normally 
printed on 5 mil thick Rogers Dtiroid 5880 dielectric substrates 
with = 2.2. Design curves have been provided, which 
appear in Figures 14-17, for structures with parameters dis- 
cussed above for the 26.5 - 40 GHz, 40 - 60 GHz, 60 - 90 GHz 
and the 90 - 140 GHz waveguide bands. 

It can be seen from Figures 14-17 that even with fin 
separations of a few mils, low values of impedance are diffi- 
cult to achieve. In some applications lower impedances than 
those obtained here which are on the order of 125 - 150 ohms, 
are desirable. One method of accomplishing this objective 
would be to increase the dielectric constant. This can be 
seen by comparing Figure 12 for a W/D = 1 and Figure 15 
for W/b = .053. For this W = 5 mil slot, it is found that 
the impedance decreases from about 165 ohms to 90 ohms when 

the dielectric constant is increased from e = 2 . 2 to e =20. 

r r 
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Figure 15. Wavelength Ratio A and Characteristic 

Impedance Z vs. Frequency for Fin-line with 
WR(19) Shield. Fins are Centered 
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Figure 17. Wavelength Ratio \'/\ and Characteristic 
Impedance Z vs. Frequency for Fin-line 
with WR(8) Shield. Fins are Centered 
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The disadvantage of this method is, of course, a reduction 
in the wavelength by a factor of 2.5 causing the already 
small circuit dimension to becom.e even smaller. 

Another possible way of achieving a lower value of impe- 
dance may be in the relocation of the fins toward the side- 
walls of the guide. Figure 18 illustrates the v/avelength 
and the characteristic impedance of a fin-line with the fins 
shifted midway between the center and the sidewall of a WR(19) 
waveguide. A comparison of Figures 18 and 15 indicates a 
2-4% change in the wavelength of the guide for W/b = 1 and 
less for smaller values of W/b. 

The values of the characteristic impedance for small 
values of W/b are relatively .unchanged since the fin-line 
structure behaves like slotline and the presence of the shield 
has little affect. For values of W/b > 0.5, noticeable 
decreases in the impedance have taken place. The most signi- 
ficant decrease occurs for W/b = 1 where the impedance 
decreases from 500 ohms for the centered dielectric slab to 
about 275 ohms for the off-centered configuration. In this 
condition the structure is a slab loaded waveguide and since 
the low dielectric constant results in little change in the 
empty guide fields, the results are as expected based on the 
voltage impedance definition for the empty rectangular wave- 
guide. A point of interest is the impedance for W/b = 1 
falling below that for W/b = 0.5. This indicates that for 
the off-centered configuration, the impedance increases with 
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Figure 18. Wavelength Ratio X*/\ and Characteristic 
Impedance Z vs. Frequency for Fin-line 
with WR(19) Shield. Fins are Located 
Half Way between Side Wall and a Guide 
Center 



W/b until a maximum is reached and then decreases until 
W/b = 1. Overall, the relocation of the fins toward one side 
of the waveguide appears not to have achieved the lowering 
of the impedance for small values of W/b as anticipated. 
Therefore, of the two methods discussed, increasing the 
dielectric constant accomplishes the lowering of the char- 
acteristic impedance desired; but as mentioned before, the 
circuit dimension shrinkage is unavoidably linked to this 
action . 
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VI. CONCLUSIONS 



A spectral domain technique for analysis of the dominant 
mode of fin-line has been presented. It has been shown that 
this efficient numerical method can be implemented compu- 
tationally to accurately calculate the dispersion character- 
istic and characteristic impedance as a function of frequency. 
A matrix formulation of the problem permitted the elements 
of the dyadic Green's function to be computed and circumvented 
the extensive algebraic manipulations associated with the 
formulation of the resulting eauations. Nximerical results 
obtained by using this method have been presented and compared 
with existing data for various structures. Excellent agree- 
ment was obtained in all cases thus establishing accuracy of 
the method and the versatility of the computer program. 

Design curves for mil lime ter -wave fin-lines of practical 
interest were also included here. Both the centered and 
off-centered fin configurations were discussed and the off- 
centered fin location was shown to result in no significant 
impedance change for small values of W/b. 

In this thesis, particular interest was devoted to the 
analysis of fin-line. The method employed, however, is com- 
pletely versatile since it is a general formulation of an 
inhomogeneous transmission media problem. Therefore, it can 
be clearly seen that by varying the parameters of the fin-line 
structure, it may exhibit the characteristics of ridged 
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waveguide, dielectric slab loaded waveguide, slotline, and 
conventional rectangular waveguide. All of these structures 
are fin-line substructures and as such are included in the 
approach presented. 
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APPENDIX A 



Derivation of Spectral Domain Matrices 
[Mg] and [Mj] 

The continuity conditions given in equations (9) through 
(18) are transformed via (24) into the two dimensional 
Fourier domain which yields 
at V = D 
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at y = 0 
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where y = = ^2 ~ >^ 3 * 

The solutions to the two Helmholtz equations given by 

equations (26a) -(26f) are now substituted; the dependence 
. 0 h 

of the coefficients A through D on a ^ is understood and 
from here on is not shown explicitly. 

(k^j^sinhYj^hj^) A®- (k^2sinhY2D) B®- (k^2COShY2D) = 0 (A-3) 

- ( j a^rsinhYj^hj^) A®+ (j a^rsinhY 2 D) B®+ ( j a^rcoshY 2 D) C®+ ( j uy Yj^sinhY^h^^) A^ 

Vi In 

+ ( jcjy Y 2 ^oshY 2 D) b“+ ( j coy Y 2 sinhY 2 D) C = 0 (A-4) 

(k^2)C^"<^c3 sinhY3h2)D® = 0 (A-5) 

- ( j a^r) C®+ ( j a^rsinhY3h2) D®- ( j coy Y2^ (A- 6 ) 

+ ( j (oy Y2sinhY2h2) D = 0 

(k^2)C^ - (k^3COshY3h2)D^ = 0 (A-7) 
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(A- 8) 



( jue2Y3COShY2h2) D®- ( ja^p) 

■U 

+ (ja rcoshYoho)D = 0 

n o z 

0 h 

- ( j a^rsinhYj^hj^) A + ( j wu Yj_sinhY^hj^) A = (A-9) 

(k^^sinhYj_h^)A® = (A-10) 

(k^lCOshYihi) A^-(k^2sinhY2D) B^- (k^2^oshY2D) (A-11) 

- ( j (jj£ j^Yj_coshYj^hj^) A®- ( j toe2Y2COshY2D) B®- ( jaje2Y2SinhY2D) C® 

(A-12) 

-(ja rcoshYi ) a“+ ( j a rsinhYo^) B + ( j a rcoshYoD)c” = J (a ) . 
n il n Z n Z zn 

This set of linear equations (A-3) through (A-12) may 
be written in matrix form, such that 
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where the 8x8 matrix is [M^] of equation (29a) , obtained 
by combining (A- 3) through (A- 10) . 

The combination of equations (A-3) through (A-8) , (A-11) 
and (A-12) form, the matrix equation (29b) repeated below in 
expanded form. 
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The matrix elements of [M^] and [Mj] are normalized at 
this point with respect to D, the dielectric substrate 
thickness; but before describing them, the normalized con- 
stants used in these elements are presented. 

(y^D)^ = (Y3D)^ = (aj^D)^+(2TT)^[(X/X')^ - 1] (D/X)^ 

(Y2D)^ = (27t)^[(X/X’)^ - Sj.] (D/X)^ 

2 

If (y^D) < 0 then the following is defined 
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= -(Y3_D)^ = -(Y3D)^ 

(Y2"D)^ = -(Y2D)^ 

Other normalized constants are 

(k ,D)^ = (k^.D)^ = (2tt)^[1-(A/X')^] (D/A)^ 
Cl c J 

(kc2D)^ = (2 tt)^[e^-(A/A')^] (D/A)^ 

omjD = 240 Ti^ (D/A) 



coe^D = cje3D = 1/60(D/A) 



toE2D = e^/60(D/A) 



SD = 2tt (D/A) (A/A ' ) 



a D = 
n 



n2TT (D/b) 



(|) even 



(2n-l) TT (D/b) cf) odd. 



The normalized matrix elements are now presented in two 

2 

forms. The first of the element equations is for (Yj^D) > 0 

2 

and the second is for (Yj_D) < 0. For the matrices [M^] 
and [Mj] , the elements m^^^^ through m^g are the same. 
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(a^^D) (6D) cos[ (y 3 _"D) (hj_/D) ] 






I 

I 




J 

"'se 



j 

m„ 



J 

^8d 



(a^D) ( SD) cosh [ (Y]_D) (hj_/D) ] 
(a^D) (3D) COS [ (yj_"D) (hj_/D) ] 
- (^j^D) ( 3D) sinhY2D 

< 

-j (a^D) ( 3D) sinY2"D 

0 



Once the matirx elements are determined for the existing 

condition (y^ < 0) , the [M^] ^ is created for the purpose of 

0 h 

eliminating the unknown coefficients A through D in (31) . 
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APPENDIX B 



Time Average Power Flow in the Spectral 
Domain of Fin-line 



The following expressions for the time average power flow 

in the three spatial regions are derived from equation (40) . 

The solutions to Helmholtz's equations given in (26a) through 

(26f) are substituted in the transformed scalar potentials 

of (41) . In the following derivations the dependence of the 
0 h 

coefficients A through D on the variable is understood. 
Also to be noted is that = Y 3 since regions 1 and 3 are 
assumed air filled. 

Region 1 

2 

Case la: (Y]_D) ^ 0, y^_ 

= A® sinh Y]^(D+h^-y) 



The term P, will be used for the power flow in region 1 for 






0 

= -Yj^A cosh Yj_(D+h^-y) 





2 

the case (y^^D) > 0. Substituting the Helmholtz equation 

solutions for region 1, (40) becomes 
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i 



D+h, 



^ ^ ^ / { [-Bu)£^a^|A®|^-6wy^Y^lA^|^+j6^ot^7 



2 i,h|2, 



la 



n=-co D 



■^^l“n^l [sinh^Yj_ (b+hj_-Y) ] 



r n 2|,h|2 - 2| S|2 .q 2 a®/7\^\* 

+ [-Bojy^a^lA 1 -Booe^Y]_|A [ +;] B (A )* 



-jkJanY;i^A^(A®) *] [cosh^Y]_ (D+h^-y) ] }dy . 



Carrying out the integration 



D+h, 

2 

/ sinh Yj_(D+h^-y)dy 



•^^[sinh (2Yj_h^) - 2Yj_h^] 



D+h, 



/ cosh^Y 3 _ (D+h^-y) dy = j^[sinh (2Y^hj^) +2Y^h^] 



Equation (B-1) becomes 



“ Bu) £, a‘ 



= -^ Re I {[ J^|A®|^-Bo)y^YilA^|^+jB^a^A®(A^)* 



la 8b 



n=-oo 



- jkJa^(A®) *A^] [sinh(2Y3^h^)-2Y3^h^] 



+ [- ^‘"^l'"ri |A^|^-3a)£j_Yi|A^|^+jB^a^A^(A^)* 



■jk^a^A^(A®) *] [sinh (2 y j_h^) +2Yj_h^] } 



;,^A®(A^) * 



(B-1) 



(B-2) 
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I 

I 

I 

I 

I 

I 

I 

I 



Equation (B- 2 ) is the final result for the power flow 

2 

in region 1 when (yj_D) ^ 0 . Putting (B- 2 ) in normalized 
forms produces the following result. 



la 8'b 



1 n ” ((^e-|D) (a D) a® 9 9 

- Re I {[ |^|^+(BD) (cjy.D) (yiD) 1 ^ 1 ^ 



n=-oo 



(YjlD) 



1 -^2 



2 * 2 (A^) * A^ 

-j (6D)^(a D)— (k D)^(a D) ^ 



h^ h, (6D) (o)U,D) (a D) 9 

[sinh2(YiD) (-A)-2 (yiD) (-^)]t[ ( T ^ d— — l^r 



A® I 2 



0 h * 

+ (BD) (o^e.D) (y,D) l^r-j (BD)^(a^D)^ 

D D D 



h s * h h 

+ j(k^D)^(a^D)^ [sinh2(Y;L°^ (-^)+2 (y^D) (-^) ] > (B-3) 



2 

Case lb: (Y]_D) < 0, yi is imaginary in which case 

(yi"D)^ = -(y^D)^ 

= jA® sin Yjl" (D+hji^-y) 

4>^ = A^ cos Yj^" (D+hj^-y) 

34 >® 

3^ = -jYjL"A® cos YjL"(D+h^-y) 
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w 



= Yi"A^ sin Yi"(D+h^-y) 



Substituting the values into equation (41) , integrating and 
normalizing with respect to D results in 



lb 8 'b 



1 ~ (SD) (cae,D) (a^D) ,e . 

- i(g) Re I {[— 1 n ,A ,2 



n=-oo 



(Yi"D) 



-^l‘-+(BD) (o)y,D) (y, "D) I^I ^ 
D 



2 (A^) * 2 (A^) * A^ 

-(6D)^(a D)% iA_; A j 

D D D D 



hj^ (SD) (tjjyj^D) (a^D) ^h 2 



[2{y^"D) (-^)-sin2(Y^"D) (-^)] + [ ^ 



0 0 h * 

A^i2, 2 , ^.A® (A^) 



+ (SD) (o)e,D) (yt ”D) | ^T+ ( gD) "• (ct^D) % 

D ^ D D 



0 * h h h 

+ (k^D) ^(a^D) ^^ 2 ^ ^][2(y^"D) (-^)+sin2(Y;L"^^ ^ ^ 



This is the final result for the power flow in region 1 
2 

when (yj_D) < 0. 

A shorter version of the derivation will be used for the 
remaining cases wherein the regional cases will be described 
and the final normalized result will be shown. This will be 
similar to that done for case lb. 
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Region 2 



For region 2, the power flow expressions are the same as 
for unshielded slowline except that the Fourier integral is 
replaced by a summation and the interval 2tt is replaced by 
the interval b. 

2 

Case 2a: ^ ^2 teal 



0 0 
B sinh YoY + C cosh YoY 



h h 

B sinh Y 2 Y C cosh y 



9 $; 

W 



0 0 
y cosh y 2 Y + C sinh Y 2 Y] 



9 $, 



h 

y ^ [B cosh y 2 Y + sinh Y 2 Y] 



Substituting into equation (41) and integrating followed by 
normalizing with respect to D produces the following. 
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I 
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II 

I 

li 

( 

i 

I 

I'' 

il 



2a 



1,D, . ; 

8*b’ A (y,D) Ip2l * (y^D) 1^2 



n=-oo 



(Y 2 D) 



9 9 

+ (6D) ii^e^D) iy^D) |!^p+(SD) (y^D) |^| 



.e ,^h 



+ j (6D) ^(a D) i^^+i^-^)-j (k D)^(a D) 



^ 2 ^2\.2 2 



D 



D D D 



[sinh(2Y2D)-2Y2D)]+[ 



( 6D) (co ^2^) c® 2 



(Y 2 D) 



D 



(6D) (coy2D) 2 



B® I 2 



(Y2°) 



y| +(6D) (cd£-D) (YoD) |-^ 
D ^ ^ D 



+ (6D) (coy.D) (y.D) j^T+j (6D) ^(a D) ) 

^ ^ D ^ D D D D 

h 0 * e * h 

-j (k_D) ^(a^D) + %] [sinh(2Y.D)+2Y.D] 

(a^D) +(Y2D) ^ ^gS .ge^* ^e 



+ [ (6D) (cj£2°) ( (y d) 



-) ( 



rY2 ^2 

D D 



2 2 
D 



ifPDW ■^^^2°^^. .B^ (C^) * (B^) * 

+ j2(6D)^(a D) (^ 11^+^ g)2(^ 



.((CV C^,(bV B^it,,3h(2Y2D)-l]> 

D D D D 



(B-5) 



This is the final result for the power flow in region 2 when 



(Y 2 D) 1 0. 
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Case 2b: 



2 

(Y 2 D) < 0, Y 2 imaginary in which case 

(y2’’D)^ = 

$2 = sin(Y2"y) cos(y2"y) 

$2 = jB^ sin(Y2"y) + cos(Y2"y) 

ai>® 

^ = jV2”B® cos yj-y - y^" sin y^-y 



= jyj"®*' cos y2"y - yj"C*' siny2"y 

2 

The final result for region 2 when (y 2 ^^^ after integration 

and normalization follows. 
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i 



P„. =-|(g) Re I {[ 



(6D) (we2D) (a^D) ^ ^ ( BD) (u)U2D) (a^D) gh ^ 



2b 



n=-oo 



(Y 2 "D) 



+ • 



(Y 2 "D) 



e h 

+ (BD) {(o£2D) (Y2"D) |%1^+(8D) {u^p^D) iy^"D) 

.(6d)2(. d)(5j 

"" 

+ (k D)^(a D) (- ^g - p- ] [2y "D-sin(2Y2"D) 1 

(BD) (u)£-D) (a^D) ^ . (BD) (uy.D) (a^D) ^ „ 

+ r 2 ^ _ ,. I I 2 ^ 2 n 1^1 2 

^ (Y2"D) 

e h 

+ (BD) (o)£2D) (Y2"D) 1^|^+(BD) {(0^2°) 

Vj*pp h* h p*.p*h 

.,6D)2(„ D, ,<B^ c!.B_! )2, (C^.rn) C , , 

^ D D D ^ ^ D D 

(a D)^-(y,"D)^ 

• [2Y2"D+sin(2Y2"D) ]+j [ (BD) (u)£2D) ( ■ ("^■■d) '> 

B® (c^) (B^)* C® U D)^- {y "D)^ 

%) + ( 8 D) (ojy^D) (-^^^ 7 — n^T ) 

D D D D ^ ^^2 ' 



hh*h*h ph*©}i 

D D D D D D D D 



0 * h 0 * h, 

+ 2(k,D)^(a^D) ® ) ] [1 -cos(2y,"D) ] } (B-6) 

^ ^ D D D^ D 
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Region 3 



Case 3 a: (Y3D) i. 0 , Y3 is real 



D sinh Y3(h2+y) 



D cosh Y3(h2+y) 



9 d>; 

W 



Y^D cosh Y3(h2+y) 



9 4 >. 

W 



Y^D sinh Y3(h2''‘y^ 



Substituting these expressions into ( 41 ) , integrating and 
normalizing yields 



3 a 



T n ” (SD) (coe^D) (a^D) ^e - 

-i(g) Re I {[— 3 n ,D ,2 



n=-oo 



(Y3D) 



^^1 "+(6D) (com^D) (yoD) 

D 



+ j(6D)^(a D) % j (k D)^(a 

h- h_ (6D) (wy-D) (a^D ) ^ . 

.[sinh2(Y3D) (-^)-2 (y3 D) (y)] + ( 



D® I 2 . . 



0 h * 

+ (gD) (o)e.D) (y.D) (gD) ^ (a„D)y 

2 3 'j^2' n j^2 ^2 



0 * h 

-j (k3D) ^(a^D)I^y- [sinh2(Y3D) (-^)+2 (y3D) (-^) ] } 



(B- 7 ) 
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This is the final result for region 3 power flow for 

(Y3D)^ ^ 0. 

2 

Case 3 b: Y3 is imaginary in which case 

(Y3"D)^ = -(y3D)^. 

$3 = jD® sin Y3"(h2+y) 

$3 = cos Y3" (^2+7) 

3^ = jY3"D® COS Y3"(h2+y) 

9 

^ = -Y3"d'’ sin Tj'Mhj+y) 



Substituting into equation ( 41 ) , integrating from y = -h2 
to y = 0 and normalizing the result with respect to D produces 
the following. 



1 n ” (u£.,D) (a^D) r,e _ 

'“ 3 b = - 5 <E> ^ ft- ' " 



n=-oo 



(Y3”D) 



i^r+(6D) (o)P,D) (y."D) 

D 



+ (6D)^(a D)^ I 2 ^+(k D)^(a 



^2 ^2 
• [2(Y3"D)(-^)-sin2(Y3"D) (-^) ] + [■ 



(6D) (o)y3D) (a^D) ^ 



(Y3"D) 



iD® , 2 



e h * 

+ (SD) (ue^D) (y"D) "■- ( 6D) ^ (a^D) ^ 



2 ^2 ^2 
-(k3D) [2 (y 3"D) (-^+sin2 ( y3"D) (-^) ] > 



D D 



(B-8) 
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This is the final normalized expression for the power flow 

2 

in region 3 under the 0 condition. It may be 

noted that a minus sign appears in front of all the equations, 
this is because a wave traveling in the -z direction was 
assumed. 
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APPENDIX C 



Computer Program 'FIN-LINE' 
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SECOND CARD: SP^CIFTPS THE NUMBER OF INITIAL VALUE DATA CAROS 

THAT FOLLOW. MAXIMUM OF 9 CARDS ALLOWED. 

COL. l: ICAR (II FORMAT) 
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IMPLICIT INTEGERd-J), REAL ( A-H, K-Z ) 
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ITERATICN LOOP FOR THE DISPERSION CHARACTERISTIC 
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IF(TEMP) 330» 720, 330 
DIVIDE PIVOT ROW BY PIVOT ELEMENT 
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